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Abstract: 

Recent measurements of the electron and muon neutrino masses squared are in- 
terpreted as an indication that neutrinos are faster than hght particles - tachyons. 
The tritium beta decay amplitude is calculated for the case of the tachyonic electron 
neutrino. Agreement of the theoretical prediction with the shape of the recently 
measured electron spectra is discussed. Amplitude for the three body decay of the 
tachyonic neutrino, Vi —>■ UiUjVj, is calculated. It is shown that this decay may 
explain the solar neutrino problem without assuming neutrino oscillations. Predic- 
tions for short and long baseline experiments are commented. Future experimental 
activities are suggested. 
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1 Introduction 



In all recent measurements the values of the mass squared obtained for the electron []T], ^, ^ ||, 
H H S Hi muon 1^, [l^ neutrinos are negative, as summarised in table 0. This observation 
requires an explanation which should be searched for on the grounds of both conventional as 
well as nonstandard concepts. 

At present the best method of determining the mass of the electron neutrino is to fit an 
appropriate theoretical formula in the endpoint region of the electron energy spectrum measured 
in the /3 decay of tritium (electron kinetic energy at endpoint amounts to ~ 18570 eV). 
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Table 1: Recent results of the electron and muon neutrino mass squared measurements. The 
first uncertainty is statistical, the second systematic. 



Improving resolution of spectrometers used in certain recent experiments lead to an obser- 
vation of two unexpected details in the electron energy spectrum, presently often referred to as 
anomalies. 

One anomaly appears as an enhancement located close to the endpoint in the integral 
electron energy spectrum [Q, 0, ^, as shown in figs. |I] and |^. Branching ratio of this structure 
comes out to be of the order of 3 ■ 10~^ or 10~^° |Q . Presently an effort is spent to understand 
the effect on the grounds of atomic physics |jll[ and statistics [|12| . 



The second anomaly, located at lower energies, is related to the slope of the measured 
electron energy spectrum. If the theoretical spectrum is fitted in the energy range below 
approximately 18300 eV then the endpoint occurs at a lower value of energy. This anomaly 
was interpreted as a manifestation of a missing spectral component with about 4% branching 
ratio and the endpoint energy around 18400 eV or 18500 eV [^f. 

The mass squared of the muon neutrino has been determined from the measured muon 
momentum in the vr"*" decay at rest. The latest, most precise results also indicate a negative 
value . The two entries in table correspond to different values of the charged pion mass 
T0| . This ambiguity soon should be resolved by additional measurements P|. 



Until now there have been no experiments dedicated to determine the tau neutrino mass. 



The latest analysis [13| lead to a result with large uncertainties and as such cannot be considered 



in the context of the idea presented in this paper. 
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Figure 1: Ratio of measured and fitted spectra (Kurie plot) from the LLNL experiment 




Figure 2: Integral electron energy spectrum in tritium decay (Kurie plot) measured in tlie 
Droitsk experiment 1^. ' 
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In view of the above results for the neutrino masses squared we elaborate on the hypothesis 
that neutrinos are tachyons. The aim of this paper is to present a confrontation of experimental 
data with theoretical calculations performed for processes involving tachyonic neutrinos. It will 
be shown that the measured shape of the electron spectrum near the endpoint in the tritium 
decay is qualitatively consistent with that predicted from the decay involving the tachyonic 
electron neutrino. The amplitude for the three body lepton flavour conserving decay of the 
tachyonic neutrino, z/j — UiUjVj, is calculated and predictions for certain phenomena discussed. 

In the following sections various values of tachyonic neutrino masses appear as a result of 
certain assumptions. Tachyons, like massive and massless particles, are subject to gravitational 
interactions too. However at the present, early stage of our work we have not yet grounds to 
comment possible constraints on tachyonic neutrino masses resulting from considerations on 
the cosmological scale. 



2 Theoretical considerations 

It has been a common conviction that tachyons could not exist for several reasons related to 
fact that they are moving with superluminal velocities: negative energies, causality violation, 
vacuum instability and others. Despite a significant theoretical effort no solution of these prob- 
lems has been found within the framework of the Einstein-Poincare (EP) relativity. However 
the recently proposed causal theory of tachyons [^4], |T^ is free of these difficulties. 



The main idea is based on two well known facts: {i) the definition of a coordinate time 
depends on the synchronisation scheme; {ii) synchronisation scheme is a convention, because 
no experimental procedure exists which makes it possible to determine the one-way velocity 
of light without use of superluminal signals. Therefore there is a freedom in the definition of 
the coordinate time. The standard choice is the Einstein-Poincare (EP) synchronisation with 
one-way light velocity isotropic and constant. This choice leads to the extremely simple form 
of the Lorentz group transformations but the EP coordinate time implies a covariant causality 
for time-like and light-like trajectories only. We choose a different synchronisation, namely that 
of Chang-Tangherlini (CT) preserving invariance of the notion of the instant-time hyperplane 



T7| , |T8[. In this synchronisation scheme the causality notion is universal and space-like trajec- 
tories are physically admissible too. The price is the more complicated form of the Lorentz 
transformations incorporating transformation rules for velocity of distinguished reference frame 
(preferred frame). The EP and CT descriptions are entirely equivalent if we restrict 
ourselves to time-like and light-like trajectories; however a consistent description 
of tachyons is possible only in the CT scheme. A very important consequence is that 
if tachyons exist then the relativity principle is broken, i.e. there exists a preferred frame of 
reference, however the Lorentz symmetry is preserved. As we know, in the real world we have 
such a locally distinguished frame - namely the cosmic background radiation frame. 

The interrelation between EP {xe) and CT (a;) coordinates reads: 

x% = x^ + u^ux , Xe = X (1) 
where is the four- velocity of the privileged frame as seen from the frame {x^). 

In refs. |15|, |T6[ a fully consistent, Poincare covariant quantum field theory of tachyons 
was given. In particular, the elementary tachyonic states are labelled by helicity. In the case 
of the fermionic tachyon with helicity —\ the corresponding field equation reads: 

(7^ (z7(9) - = (2) 
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where the bispinor field ip is simultaneously an eigenvector of the helicity operator with eigen- 
value Here the 7- matrices are expressed by the standard ones in analogy to the relation (|Tp. 
The solution of (||) is given by: 

^(x, u) = f (fk 5{e + K^)9{k'>) \w{k, u)e"'%\k) + v{k, u)e'"'''a{k)] (3) 

(27r)2 J L J 

where the operators a and b correspond to neutrino and antineutrino, respectively. The ampli- 
tudes V and w satisfy the following conditions: 



w{k, u)w{k, u) = {k- 7^/c7)^ - ^J^^^^j, 
v{k, u)v{k, u) = -{k + f A;7)^ (1 - 1^7=== 



(4) 



(5) 

w{k,u)'y^U'jw{k,u) = v{k,u)'y^U'jv{k,u) = 2q (6) 

w{k^, u)'y^u-fv{k, u) = 0. (7) 

Here q = u^j^k^^ is equal to the energy of the tachyon in the preferred frame and H denotes 
the space inversion operation. It is easy to check that in the massless limit k — >^ the above 
relations give the Weyl's theory. 

3 Tritium (3 decay with the tachyonic neutrino 

Until now the mass squared of the electron neutrino has been determined by fitting the mea- 
sured integral electron energy spectrum with a function comprising an expression for the decay 
probability. The integral form of this expression is as follows (neglecting summation over final 
states) 0: 



S{E) ~ pe(^e + me)(^o - E,)J{E^ - - m- 



where i?e, Pe, f^e are the electron kinetic energy, momentum and mass, Eq is the endpoint 
kinetic energy and m^^ is the square of the electron antineutrino mass. This formula is obtained 
from the standard calculation of the three body weak decay, including a massive neutrino. 
However if one wants to test the hypothesis that neutrinos are tachyons, this formula is no 
longer valid for that purpose even if the sign in front of the m^^ term is changed. The correct 
expression can be obtained from a similar, although by far more complicated, calculation. 

Let us consider the j3 decay in the framework of section 0, using effective four-fermion 
interaction. In the first order of the perturbation series the decay rate for the process: n — > 
+ e~ + Ve with the tachyonic electron antineutrino reads: 
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where (i$3 is the phase-space volume element: 



while 



d$3 = ^(A;O)^(/O)0(r°) 

6{e - ml) 6{l^ - ml) 6{r^ + k^) 6\p - k - I - r) d^k dH d^r 

= 2 Tr [UeUel^WW'^''] Tr UpUp'^^ (l - QaI^) UnUnlv (l - QAl^ 



(10) 



fill 



Here p, fc, /, r are the four-momenta of n, e~ and 17 respectively; the corresponding masses are 
denoted by m„, m^, me and k (in the case of tritium decay m„ and denote the masses of 
and ^He). Gp and (yf^i are the Fermi constant and the axial coupling constant. The amplitudes 
satisfy usualQ polarisation relations: m„m„ = P7 + m„, WpWp = k'j + mp, 



llp^ Up^ Ue 



Ij + me whereas ww is given by eq. ^j. After elementary calculations eq. |Tl] reads: 
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Because the Solar System is almost at rest relatively to the cosmic background radiation 
frame0 we can calculate the electron energy spectrum dT/df in this frame with sufficient 
precision. The resulting formula is rather complicated (|16|); we present it in the Appendix. 

Formula (^ gives identical results in the limit ^ as that for the massless neutrino. 
Differential electron energy spectra near the endpoint for massless and tachyonic neutrinos are 
shown in fig. ^ (we used the following values for the masses of and ^He: 2809.94 MeV 
and 2809.41 MeV). The distinctive feature of the spectrum near the endpoint in the tachyonic 
case is the enhancement with a step-like termination. Function (|1^) goes to zero over a narrow 
range of energy, contrary to a monotonic character in the cases of the massless (and massive) 
neutrino. In order to allow for a qualitative comparison with the measurement of ref. |Q, the 
electron spectrum given by (|l^) was integrated and smeared with constant experimental energy 
resolution of 4 eV (assumed Gaussian). 

As a result the step-like endpoint structure of the differential spectrum appeared as an 
enhancement in the integral electron energy distribution, as shown in fig. ^3. Position, mag- 
nitude and width of the enhancement depend on the tachyonic electron neutrino mass, k. At 
this stage of investigations we may conclude that using formula (0) it is possible to reproduce 
qualitatively the enhancement in the electron energy spectrum near the endpoint if a sufficiently 

"'^but with 7 in the CT synchronisation 
^relative velocity is about 370kni/s 
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Figure 3: (a) Theoretical differential energy spectrum dT/dE for tritium decay with the tachy- 
onic electron antineutrino of mass n; (b) smeared integral spectrum (Kurie plot) - endpoint 
energy indicated by arrow. Curves drawn for masses: k = 0, 10, 20 eV. 



large value of the tachyonic neutrino mass is taken. Indeed, no ''bump' is seen for k = 10 eV 
while it is already clearly marked for k = 20 eV. It is indispensable to perform a fit to the 
experimental spectrum; however we were not in possession neither of the data nor the details 
of the apparatus. 

For the purpose of considerations presented in section ^ we need an input value for the 
tachyonic electron antineutrino mass. The value k ^ 5 eV, quoted in refs. [@, @], is not 
meaningful for the tachyonic hypothesis since it was obtained by fitting function (§) which 
describes the massive neutrino {m? > 0). As an 'educated guess' (fig. §3) we take 20 eV for 
the mass of the tachyonic electron neutrino. 

4 Decays of tachyonic neutrinos 
4.1 Introduction 

Tachyonic neutrinos are in general unstable. Similarly to a three body decay of a massive 
particle, a tachyon may decay into three tachyons. In the latter case however a specific channel 
is allowed with the initial tachyon appearing also in the final state. In the case of tachyonic 
neutrinos the corresponding process would be the decay of a neutrino of flavour 'i' into itself 
and a neutrino - antineutrino pair of flavour 'j': z/j ViVjVj, where i,j = e,fi,T, as shown in 
fig- i 

Another important mode is the radiative decay, Ui Ui'-f, shown in the same figure. In 
both processes the lepton family number is conserved, contrary to those involving a massive 
neutrino. In this section we deal with the three body decay and calculate partial widths Fjj 
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Figure 4: Lepton flavour conserving decays of the tachyonic neutrino. 



in order to find out whether the tachyonic hypothesis is not in contradiction with estabhshed 
observations. 

Tachyonic neutrino decays mean an additional, quahtatively similar process to oscillations. 
Consider a beam of neutrinos of a given flavour, characterised by a certain energy spectrum. 
As a result of decays new flavours appear as a function of distance (including antineutrinos 
in neutrino beams and vice versa) but in a pattern much more complicated than in the case 
of oscillations. /^From the theoretical point of view the decay has a clear advantage over the 
latter since its magnitude can be calculated within the theory of weak interactions; masses of 
tachyonic neutrinos, which can be determined from independent measurements, are the only 
input parameters in that case. 



4.2 Derivation of formulae 

Assuming universality of weak interactions we can calculate the decay rate for the process 
i/j ^i^j^j ihj = e, /U, r) which is given by: 

dTi. = , \ \T? , (13) 

where 

rf$ = e(/)^(/0)^(r°)5(/ + fi:2)^(;2 ^ ^2)^(^2 ^ ^2)^(^ - p - / - r)d^pdHd^r (14) 

is the phase space volume element, the corresponding four-momenta k,p, I, r as shown in fig. ^ 
K and jj, are the masses of the i-th and j-th neutrino, q - energy of the i-th neutrino and 

|T|' = 4G|Tr [v,ik)v,ik)^^v,ipMp)r] Tr [w,ir)w,ir)^^v,{l)v,{l)^,] (15) 

can be calculated using results of section ^ The final formula for r(g, k, /i) is rather long and 
complicated and for this reason we do not present it here. 

4.3 Results and discussion 

Existence of the three body decay raises questions about stability of tachyonic neutrinos in 
different environments, ranging from laboratory to the Universe. As shown above, the mean 



8 



lifetime of the tachyonic neutrino of a given flavour can be calculated only if masses of all three 
neutrino species are known. The present situation in this respect may be viewed as such that 
the mass of the tachyonic electron neutrino may be considered to be approximately known 
(k ~ 20 eV) but not those of the muon and tau neutrinos. Uncartain measurements in the 
two latter cases are the principal difficulty in discussing the consequences of neutrino decays 
in detail. For that reason the following remarks are to some extent of a speculative nature. In 
order to simplify the discussion, we limit our considerations to a two flavour scenario with the 
electron, u^,, and the heavy, u^, neutrinos, implicitly assuming ^ i^h- 

It is a specific feature of tachyonic neutrinos that mean lifetime for the three body decay 
decreases as a power of energy, q. Magnitude of this dependence varies significantly over a wide 
range of tachyonic neutrino masses, as can be seen in fig. ^. In general the process of neutrino 
energy degradation is slowing down with subsequent decays. 

Partial width for the decay of a tachyonic neutrino into a given final state increases with 
the mass of the neutrino from the pair, /ij. For a light neutrino like this dependence is 
very strong in a wide range of /ij masses, as can be seen from fig. |^a. Total width for the 
electron neutrino decay {i = e, j = e,H), neglecting the radiative decay, is approximately 
equal to the sum of partial widths for the two channels: = Fee + ^eH- If i^h i^e then 
Fee ^ ^eH SO to a good accuracy Fg ~ Fg/f which means that the mean lifetime of the tachyonic 
electron neutrino is determined by the (unknown) mass of the heavy neutrino. Concerning the 
latter, the present measurements of the muon neutrino mass squared justify considering the 
value kh ~ 0.2 MeV which lies within a Icr uncertainty of the latest measurement. A heavy 
neutrino of this mass would be significantly less stable than the electron neutrino, as can be 
seen from fig. ^d; moreover, its mean lifetime does not depend on the exact value of the electron 
neutrino mass. The ratio of partial widths for both decay channels, THe/^HH, amounts to 0.1 
for g = 100 MeV, making the uh — ^ VH^eT^e decay channel not negligible. 

In general, inverse partial widths may be arbitrarily large provided tachyonic neutrino 
masses are sufficiently small. If the mass of the heavy neutrino were smaller than about 100 eV 
then the mean lifetime of the electron neutrino would exceed the age of the Universe for en- 
ergies up to those presently reached in accelerator beams. On the other hand it is interesting 
to speculate about consequences of larger masses. If the mass of the heavy neutrino were of 
the order of 1 MeV then the mean lifetime of g ~ 10 MeV solar electron neutrinos would be 
comparable to the time of fiight from the Sun to Earth (fig. |^a). Under this condition one would 
expect both fiavour composition and energy spectrum in the solar fiux to be modified due to 
decays. At a first glance this is an appealing possibility to explain the solar neutrino problem in 
a simple way, without involving the phenomenon of neutrino oscillations. However if neutrinos 
of similar (or higher) energies travel a larger distance, like e.g. from SN1987A to Earth, their 
energy spectrum would be modified too. We have not yet carried out a detailed co-analysis 
of the data from both sources. Under certain assumptions we may put an upper limit on the 
heavy neutrino mass using the SN1987A data. For that purpose we assume tentatively that the 
energy spectrum of the SN1987A neutrinos, measured on Earth, is little or not at all distorted 
as compared to the one during emission. This may be achieved if: {€) either a negligible fraction 
of neutrinos decayed during to ~ 5 ■ 10^^ s (150 000 years) on the way to Earth (if r <^ to) or 
{ii) after each decay the VjVj pair acquired a negligible fraction of the initial neutrino energy (if 
r < to)- Inequality t <^to implies the mass of the heavy tachyonic neutrino to be much smaller 
than about 5 keV, as can be read from fig. ^ for q = 0.1 -h 10 MeV. This result should be 
taken with caution since there is no firm experimental evidence for the underlying assumption. 

Three body decays of tachyonic neutrinos may simulate oscillations in terms of experimental 
signature in an experiment located at a fixed distance from the source. In order to distinguish 
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Figure 5: Dependence of the inverse partial decay width (seconds) on (prehminary results 
with accuracy of few percent due to numerical complexity of calculations). Curves correspond 
to different decaying neutrino energies, q: 0.001, 0.1, 1, 10 and 100 MeV, in the indicated order; 
(a) for the electron neutrino, k = 20 eV; (b) for the heavy neutrino, n = 0.2 MeV. 



between oscillations and decays measurements should be performed at different distances. Can- 
didate events for neutrino oscillation have been recently reported in the Vf. channel at 
LAMPF; if interpreted as such, this measurement indicates oscillation probability of 0.3% (over 
a distance of ~ 30 m) ^ . In fig. ^ we show the inverse partial decay width for the decay 



vh ^Hi^ePe cLS a fuuctiou of the decaying neutrino mass, kh-, for energy g = 50 MeV. If we 
identify with uh then it follows from the experimental conditions that l/F for this decay 
channel should be of the order of 10'^ s to observe the above fraction of decays. Such a low 
value can only be achieved for the muon neutrino mass in the range of dozens MeV, which is 
excluded by measurements. Thus the LAMPF candidate events cannot be explained by the 
three body decay of the tachyonic muon neutrino: — > z/^z/gl/g (as can be seen in the same 
figure, the result does not depend on the exact value of the electron neutrino mass). 

I 

5 Comments 

In this section we list a number of remarks concerning experimental aspects of the hypothesis 
of tachyonic neutrinos. The first and the easiest task is to reanalyse the existing tritium decay 



data by fitting eq. |T6| for the electron energy spectrum. A positive result of this procedure 
would not be conclusive for the question whether neutrinos are tachyons since one can never 
exclude a possible explanation based on conventional physics. An agreement should be reached 
concerning the shape of the experimental electron energy spectrum near the endpoint; results of 
experiments differ as to position and width of the enhancement. It is desirable to measure the 
electron energy spectrum near the endpoint and search for possible ' anomalies' in (3 decay of a 
nucleus other than tritium. One such experiment, using ^'^^Re, is in progress (endpoint energy 
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Mass of decaying neutrino, (IvIeV) 



Figure 6: Dependence of the inverse partial decay width on the mass of decaying neutrino, kh, 
for the decay into electron neutrino - antineutrino pair (fi^ = 10, 20 eV). 

m 

2.48 ± 0.04 keV) |2^. Redetermination of the muon neutrino mass squared with increased 
precision would be of great importance but very difficult experimentally. 

The theory of tachyonic neutrino decays and first results presented above open a new field 
which ought to be investigated in detail. Search for decays of hypothetical tachyonic neutrinos 
is certainly an experimental challenge. Most suited are long baseline exepriments which are 
planned in future with a purpose of searching for neutrino oscillations. Decay probability 
increases as a power of the beam energy so the highest beam energies are most desirable. The 
results presented above indicate that decays should be searched for the heavy (muon?) neutrino 
in the channel: ^^i^iJ^ii. (i-e. with experimental signature u^) which has a larger 

branching ratio than the one with the uj^f, pair in the final state. 



6 Summary and conclusions 

According to the results of several recent high precision experiments, both electron and muon 
neutrinos are found to have negative values of their masses squared. The significance of this 
statement in the case of the electron neutrino is supported by consistent results obtained in 
independent measurements. The result for the muon neutrino is not conclusive. Measurements 
of the Mainz, Troitsk and LLNL groups [l|, 0, Q reveal an enhancement in the integral electron 
spectrum near the endpoint. We calculated the amplitude for the beta decay of tritium involving 
the tachyonic neutrino in the framework of the causal theory of tachyons [p^ p!5| , p!6[] . As a 
result we reproduce qualitatively the enhancement observed in experimental spectra. 
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We presented results of preliminary studies concerning three body tachyonic neutrino decays. 
In general the decay probability increases with neutrino energy according to a power law. 
Tachyonic electron neutrinos are stable (mean lifetime exceeds the age of the Universe) in a 
wide range of initial neutrino energies. For the purpose of making reliable predictions more 
precise measurements of the muon and tau neutrino masses are needed. More studies are in 
progress. 
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7 Appendix 



The (differential) energy spectrum of electrons in (3 decay, dT/dP, may be obtained by means 
of formulae , ([ToD , ([TTD after elementary integration of: 



dr 



with 



dP 4(27r) TTl-fi J max (r-_,0) 

-Am'^f + Am'^mn + 2{Pfmn - 2fml 



dr^ 



(16) 



2(mg — 2/°m„ + m'^] 



+ 



± 



((/o)2 _ m2)((Am2)2 + 4/€2m2 - AAmH^nin - SK^^mn + An^ml + 4{Pyml 



2(ml — 2Pmn + mf^) 

where Am^ = — + ml — and the following expression for the matrix element squared: 




M(Z°,r°^ 
+ 
+ 



32 m„ TT^ 



I 2(1 - g\)mem.pK - (1 + g\)memnK 



+ 



(1 + g\)meK - 2gAK'^ - gA^m^ 
(1 + g\)m^K - 2gj,ml + ^.^Am^ 
2gAmn{ry + 2(7^m„(/°)2 + 



+ 



r0)2 + k2 



(1 + g^j^t^ ( -Am^ - mg ) + gAra^K{2¥? + Am^^ 



+ 



(1 + gA)mnK^ - (1 - Qa)^?!^^ - 2gAmemnK 



p + 



+ 
+ 
+ 



-(1 - g\)mpAm'^ - (1 + 9\)mnfi'^ - 2gAmemn 



-(1 + g\)Am'^ + 2gAmeK - (1 - g\)mn'mp 



K 
0^0 



fr°)2 + 



2^^me«: - (1 - ^7i)m„mp + (1 + ^i) (^^Am^ - 
(1 + ^i)«:^(/°)^ - (1 + gDm^Air - (1 + ^?i)m„(r°)^}} . 
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